Modeling-Error-Driven Performance-Seeking Direct 

Adaptive Control 


Nilesh V. Kulkarni 1 

Perot Systems, Inc., NASA Ames Research Center, Moffett Field, CA - 94035 

John Kaneshige 2 and Kalmanje Krishnakumar 3 
NASA Ames Research Center, Moffett Field, CA — 94035 

John Burken 4 

NASA Dryden Flight Research Center, Edwards, CA - 93523 


This paper presents a stable discrete-time adaptive law that targets modeling errors in a 
direct adaptive control framework. The update law was developed in our previous work for 
the adaptive disturbance rejection application. The approach is based on the philosophy that 
without modeling errors, the original control design has been tuned to achieve the desired 
performance. The adaptive control should, therefore, work towards getting this performance 
even in the face of modeling uncertainties/errors. In this work, the baseline controller uses 
dynamic inversion with proportional-integral augmentation. Dynamic inversion is carried 
out using the assumed system model. On-line adaptation of this control law is achieved by 
providing a parameterized augmentation signal to the dynamic inversion block. The 
parameters of this augmentation signal are updated to achieve the nominal desired error 
dynamics. Contrary to the typical Lyapunov-based adaptive approaches that guarantee 
only stability, the current approach investigates conditions for stability as well as 
performance. A high-fidelity F-15 model is used to illustrate the overall approach. 


introduction 

^ table direct adaptive flight control using neural networks received significant attention after the work of Rysdyk 
k-Jand Calise 1 . The Rysdyk and Calise baseline control design consists of dynamic inversion with command 
augmentation. The augmented command signal is given by an adaptive neural network, whose weights are updated 
online. The update law guarantees that a Lyapunov function, which includes the system error and the neural network 
weight error, is non-increasing with time. Prior to this work, various researchers, including Narendra 2 ' 4 , have used 
adaptive augmentation of the command signal in a dynamic-inversion or feedback-linearization-based control 
design. The subsequent work has looked at different aspects of the problem including stability with different kinds 
of networks, actuator saturation, and output feedback. These developments have also been implemented in various 
flight applications. 

Complementing this effort, researchers at NASA Ames Research Center have been looking at different adaptive 
control approaches within the Intelligent Flight Control program 5 ' 8 . Several approaches have been considered in the 
various guidance and control loops. These include the work of Rysdyk and Calise for inner-loop control and 
adaptive critics 8 for outer-loop guidance. While guaranteeing stability, achieving good performance has always been 
one of the major concerns in the inner-loop architecture. The key issue in all the adaptive control literature is the 
definition of the Lyapunov function. Once the Lyapunov function is prescribed, the update law guarantees negative 
semi-definite time derivative of this Lyapunov function. Typically the Lyapunov function includes the tracking error 
and the parameter error terms. In this work, we define the Lyapunov function that captures the modeling error 
between the assumed plant and the actual plant. Guarantees of a monotonical decrease of this Lyapunov function 
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then reduce the modeling error, and thereby provide performance that would have been achieved from the original 
control system with exact knowledge of the plant. The critical point is that this is achieved in a direct adaptive 
framework rather than the typical indirect adaptive control framework. In our previous work, we investigated this 
approach for the problem of adaptive disturbance rejection 9 

With this introduction, the remaining portion of the paper is organized as follows. Section II outlines the overall 
control architecture, which includes the baseline control design and motivates the adaptive control design with a 
SISO example. Section III looks at a general state space problem, parameterization of the adaptive controller, and 
presents a normalized gradient update law. Section IV presents the results of this adaptive control implementation 
for the rate-command attitude-hold (RCAH) of an F-15 model. Finally, section V concludes the paper, and outlines 
directions for future work. 

Il.Control Architecture 

Figure 1 presents the overall adaptive control architecture. We first outline the baseline controller (without 
adaptation) to motivate the underlying control principle. This is followed by discussion and equations with the 
adaptive element. 


y„om(* + i) 



Figure 1: Adaptive Control Architecture 

The control system is given a command Tcom(^ + l) (eg. Pitch rate command from the pilot’s stick). The time index 
(k + 1) refers to the desired value at the next time index. Given the knowledge of how fast or slow the plant can 
handle such a command, it is typically taken through a second order reference model with appropriate damping and 
natural frequency to get the corresponding achievable reference command T re f(^ + 1). The controller is design to 
achieve a prescribed second order error dynamics with respect to this reference command. Let this error dynamics, 
in the discrete form, be given as: 


e(A: + l) + K Pe e(A:) + K Ie e I (A:) = 0 ) ( i) 

where e (^') = k(^) — T ref (k ) ; an( j ej(£) represents the integrated error until time index k. K Pe and K[e are chosen 
appropriately to have the desired transient response characteristics. Equation (1), with the definition of the error, is 
used to compute the control input to achieve this desired error dynamics as follows. Equation (1) corresponds to: 
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( 2 ) 


[y(k + 1) - y ref (k + 1)] + K Pe [ y(k) - y ief (A:)] + K Ic e,(A:) = 0 

u 

y(k + 1) = y iei (k + 1) + Kp e [y rel ( k ) - y(k)\ - K Ic e t (A) 

The plant output ) ! (k + 1) needs to satisfy Eq. (2) to achieve the prescribed second order error dynamics. The right 
hand side of Eq. (2) can thus be labeled as + 1) , the desired plant output. Let the plant dynamics be given as: 

y(k + 1) = / [ y(k), y{k-\),..., y(k - p v ), u(k - 1), u(k - 2), . . . , u(k - p u )] 

(3) 

+g [_y{k\ y(k - 1 y(k- p y ), u(k-l),u(k-2),...,u(k- p u )]u(k) 

We can thus invert the dynamics represented by Eq. (3) to compute the control u(k) [ 0 achieve the desired error 
dynamics [Eq.(l)]as: 


1 | Trer ( A + 1) + K ,, e [ T ref ( A) - y( A )] - K (A) j 

g [ y(k), . . . , y(k - p y ),u(k-\),...,u(k- p u )] [-/ \_y(k), y(k - 1 y(k- p y ), u(k - 1), u(k - 2), . . . , u{k - p u )] j 

(4) 

This control input, with exact knowledge of the plant (/ and g), will help achieve the desired second order error 
dynamics. With modeling uncertainties/ errors, we will not know/ and g exactly, but only their estimates given by 

the model, / and g . The adaptive augmentation is now designed to offset these modeling errors, so that we can get 
the same error dynamics or the desired performance. With the adaptive augmentation, as shown in Figure 1, the 
desired output, y d es( k + 1) , is given as: 


Tdes ( k + !) = Tref ( k + 1) + K Pe [ J ref (k) - y(k)\ - K le e! (k) - y ad (k + 1) 


(5) 


The control input is given as: 

l | Wef (* + 1) + K Pe [>w (k) - y{k)\ - K Ic e, (k) -yJk + \) 1 

g [y(k), ...,y(k- p y ), u{k- 1), . . . , u{k -p u )\ [-/ [ y(k), y(k - 1 y(k-p y ),u(k- 1), u(k- 2), . . . , u{k -p„)]| 

(6) 

To analyze the effect of this control input, we look at the modeling error, which is defined as the difference between 
the actual plant output and that predicted by the model: 

e(k + 1) = y (k + 1) - / [y(k\ y(k - 1), . . . , y(k - p y ), u(k - 1), u(k - 2), . . . , u(k - p u )] 

- g \_y(k ), . . . , y(k - p y ), u(k - 1), . . . , u(k - p u )] u(k) 

Substituting the expression for the control input, given by Eq. (6), in Eq. (7) gives: 

£ (k + 1) = y( k + 1) - { JV {k + 1) + Kp e [ y re{ (k) - y(k)} - K Ie e, (k) - y ad (k + 1)} 
s{k + 1) - y ad (A: + l) = yi y k + \)~ j ref (k + 1) + K Pe [y(k) - y K( (A:)] + \L le e v (k) 

In terms of the definition of the tracking error, Eq. (8) can written as: 
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e{k + 1) + K Pe e(k) + K Ie e! (k) = s(k + 1) - y^k + 1) 


(9) 


Equation (9) represents the critical equation of this approach. The left hand side of the equation is the desired 
second order error dynamics. The right hand side of the equation is the difference between the modeling error and 
adaptive augmentation signal input. The equation, thus, says that if the adaptive augmentation signal can learn the 
modeling error and cancel it, the error dynamics of this control loop will be restored to its desired nature. In other 
words, we will recapture the performance desired from this control loop. We, therefore, define the left hand side of 
Eq. (9) as the performance error, E. 


£(& + l) = e[k + l) + K Pe e(£) + K Ie ej(£) 

We can now form a Lyapunov function that is a function of the performance error as: 

L(k) = yE(kf 


( 10 ) 


( 11 ) 


An update law now can now be devised for the adaptive augmentation input, Tad , that forces a decreasing/non- 
increasing nature to the this Lyapunov function. 

IILParameterization and Update Laws for the Adaptive Augmentation 

Section II provided a SISO example for motivating the proposed adaptive control approach. In this section, we 
investigate two main questions. The first relates to the parameterization of the modeling error, and the second relates 
to the choices for designing stable update laws. 

1. Parameterization for a linear system: 

Consider a linear system of the form: 


x(k + 1) = A x(k) + Bu(£) 

In a manner similar to that illustrated by Eqs. (4-6), the control input is computed as: 
u(£) = B >, (£ + l) + K Pc e(£) + K Ie e, (k) - x ad (k + 1) - Ax(k)j , 


( 12 ) 


(13) 


where A and B are estimates of the system A and B matrices. If the system matrices (A, B) are known, 
there is no need for the adaptive augmentation, and the control input is computed as: 

u(£) = B ' [x ref ( k + 1) + K Pc e(A;) + K Ic e, (k) - Ax(fc)] ( 14 ) 

If these control inputs need to provide the same desired error dynamics, they need to be equated, which gives 
the form of the idealized value of the augmentation signal. 

B" 1 [ x re f( A: + l) + K Pe e ( A: ) + K ie e I ( A: )- x Id( A: + l)-Ax(^)] = B“ 1 [x ref (A: + l) + K Pe e(A:) + K Ie e I (A:)-Ax(A:)] (15) 


< d (* + !) = (i -BB - 1 )[x ref (k + 1 ) + K Pc e(£) + K Ic e#)] + (bB-'A - A)x(*) 


fl-BB ) (BB A-A) 

x ref (* + !) + K Pc e(£) + K Ic e, (k) 

A / v / J 

x(k) 
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2. Parameterization for a non-linear system affine in control: 

Consider a non-linear system that is affine in control, and whose dynamics can be written as linear in 
parameters. 


x(£ + l) = W,p,(A:) + Bu(A:) 


(17) 


where "/ is the linear dynamic weight matrix, and P / correspond to the nonlinear functions of the system state. 
The control input is computed in a similar manner as: 


u(£) = B ' x ref ( k + 1) + K Pc e(£) + K Ie e, (k) - (k + 1) - W / p / (k) ; 


(18) 


where Wy an d B are the corresponding estimates of the system matrices. By carrying out the analysis similar to 
the linear system case, the ideal augmentation signal can be computed to be: 


; d (£ + l) = [(l-BB-') (BB-'Wy-Wy) 


x ref o k + 1) + K Pc e(£) + K Ie e, (k) 

P f (k) 


(19) 


Equations (16) and (19) imply that the ideal augmentation signal can be written as: 


<d(£+i)=w;/pw, ( 20) 

with the ideal weights, W ad ; and the basis functions, P , as given in Eqs. (16) and (19). It is interesting to note that 
these are the same basis functions used in Rysdyk and Calise 1 . Thus, we can parameterize a neural network in this 
form, and compute the ideal weights iteratively using an appropriate update algorithm. 

A.Update Laws for the Adaptive Augmentation: 

Having looked at the question of parameterization, we now consider the question of a stable update law for the 
parameters, W ad Parameterizing, the adaptive augmentation signal in the form given by Eq. (20), and using the 
definition of the performance error as given in Eq. (10), Eq. (9) can be rewritten as: 

E W = £-y a c#) (2i) 

Compared to Eq. (9), this is written for a vector, E, corresponding to the general case of multiple control loops. 
Written in this form, the equation savs that w p rf ' trying to estimate the vector modeling error, £ (for all loops) with 
the adaptive augmentation signal, y a d(^) . E (£) is the corresponding error in the estimate. This error dynamics for 
the performance error, E, corresponds to a system identification like problem. This opens up a host of approaches 
for doing this online system identification. In this work, we consider the normalized gradient update approach. 

3. Normalized Gradient Update' 1 : 

Let E,(£ + l) correspond to the ; lh element of the vector performance error E(& + 1) 

"ad, represent the ; lh 

column vector of the weight matrix W ad „ which corresponds to the ideal weights that minimize the performance 
error E to A* . Similarly, let "ad, represent the ; lh column vector of the current estimate of the ideal weight matrix. 
The update law for each of these column vectors of the weight matrix is given as: 
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( 22 ) 


W ad . (& + !) = W adj (k) + 


a(k)*Ei(k)* P(fc) 
'l + p r (£)P(£)] 


a (k) corresponds to the learning rate that needs to satisfy the condition: 

0 < a(k) < 2 (23) 

If the system ex^priprir-pQ enfificie"* ipp'-^+Pit °'"'*tation, reference Till proves that this weight update law [Eqs.(22- 
23)] guarantees E t (k) — > 8* and W ad . ( k ) — > W adj Vi = 1.../ as k —> °° . Appendix A provides the stability proof for 
this update law. 


B.What happened to tracking error? 

The final part of this analysis corresponds to investigating the behavior of the system error e (^) . This work is 
motivated towards providing an update only when we see modeling error, and not just if there is tracking error. 
However, having noting that, tracking error is what is ultimately of importance. It is therefore important to analyze 
the asymptotic behavior of the tracking error given the behavior of the performance error. For simplicity, in this 
analysis we consider the case where the desired error dynamics is first order given as: 

E, (k) = e,(k + 1) - K pe e t (k) = 0 ( 30 ) 

Let < & after time k, where & is some small positive scalar. This implies 

| e , (^ + 1) - K pe _ e, (£)| < £ ( 31 ) 

From Cauchy- Schwarz inequality, 

\e, {k + 1) - K pej e t {k)\ > \e t (k + 1)| - \K pei | \e, (A:)| ( 32 ) 


Equations (31) and (32) imply: 

\ e i (k + 1)| - \K pei | |e f (k)\<8 
\ e i(k + 1)| < \K pe . | |e,-(&)| + 5 
\ e i{k + 2)\<\K pe \y i {k + \)\ + S 
\e i (k + 2)\<\K pei \(\K pei \\e i (k)\ + s) + S 
<\K pe $\eAk)\ + \K pe \8 + 8 
\e, (k + 3)| < \K n f \e t (*)| + 1 K pe> fs + \K„\s + S 

\e, (k + n) | < | K pei f \e,{k)\ + (l + 1 K n \ + \kJ +... + \k„ f ) 

Since ' for stable error dynamics, as k — > oo ; |e,(£)| j s bounded above as: 
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(34) 


hw|< 
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Thus, if the performance error is bounded, Eq. (34) establishes bounds on the tracking errors. A similar analysis can 
be carried out for second order error dynamics. The result summarizes that as long as the desired error dynamics 
(first or second order) is stable, the tracking error will be bounded above given bounds on the performance error. 

IV.Application to Aircraft Control 

The modeling error-driven performance-seeking adaptive control design was implemented for aircraft roll, pitch, 
and yaw rate control. The NASA Intelligent Flight Controller (IFC) design has been tested, and is currently 
undergoing various modifications for being flight-tested on the research F-15 aircraft. The IFC design has been 
implementing the adaptive control design as outlined by Calise and Rysdyk. For implementing the performance- 
seeking adaptive augmentation, the requirement was that it needed to fit within the existing architecture. The main 
issue in the implementation is that the baseline controller in the IFC architecture uses continuous -time aircraft 
dynamic inversion, whereas the proposed design has been outlined in the discrete-time. The equations outlined in 
sections 11 and 111 have been formulated for a discrete-time model inversion. We realized, however, that after 
reducing the problem to the core error dynamics, the problems became identical. 

The error equation for the continuous-time implementation for a desired second-order error dynamics is given as: 

e + K p e + K I \edt = £-U dii (35) 


The error is defined in the same manner as the discrete case (eg. ^ 7 ,-< 7 ). The modeling error, e , however 

corresponds to the difference in the acceleration as predicted by the model and that actually observed. Similarly, 

^ad represents the augmentation acceleration command given by the adaptive block. If the left hand side of Eq. (35) 
is discretized while maintaining the continuous-time constants, the resulting discrete -time equation is given as: 

~{ e (£) + (K p \t - \)e(k - 1) + KjMe^k- 1)} = [s(k) - U aA {k)\ ( 36 ) 

Defining the left hand side of Eq. (36) as the modified performance error, E(k ) : we get 

E(k) = s{k) - C/ ad (k) (37) 

This modified performance error equation is identical to the discrete-time version given by Eq. (21). The adaptive 
augmentation acceleration signal, ^ad (^) , can be parameterized in a similar manner, and the same update laws 

remain valid for the parameters of this augmentation signal for reducing E(k ) . A zero value of this modified 
performance error restores the second order error dynamics (LHS of Eq. (35)) to zero, and thereby regains the 
desired performance from the control loops. 

Formulated in this manner, this adaptive approach fits within the existing IFC framework, and is considered as 
an alternate approach for flight testing. In the following discussion, we present the results of this implementation on 
the high fidelity model of the modified F-15 aircraft used at the Dryden Flight Research Center. In this study, we 
look at 2 cases. In the first case, the right stabilator is locked at 4 degrees at 10 seconds into the flight experiment. In 
the second case, the canard multiplier is set at -1, again at 10 seconds into the flight experiment. The behavior of the 
aircraft and update algorithm is examined for the the longitudinal and lateral pilot stick inputs given by Figure 2. 

Figures 3-7 present the behavior of the aircraft and the neural net signal for the right stabilator failure case. In 
figure 3, we can observe the learning in the pitch channel, when the aircraft pitch rate starts following the desired 
pitch rate. Figures 4 and 5 examine the learning in the roll and yaw channels respectively. The learning in the yaw 
channel is not as good as in the pitch and roll channels. Figure 6 examines the behavior of the performance errors in 
each of the three axes as a function of learning. These performance errors drive the updates in each of the three axes. 
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We can note that these get smaller with time. Figure 7 presents the aircraft surface commands for this maneuver. 

In a similar manner, figures 8-11 present the results for the incorrect canard multiplier case. Again we can note 
the good learning achieved in the pitch and roll channels. The insufficient learning in the yaw channel is being 
examined based on its baseline control system. 


Pilot Stick Input for Mach 0.75, Alt 20K - Right Stab Lock at 4 deg - Failure Inserted at t = 10 sec 



time (seconds) 



Figure 2: Pilot Longitudinal and Lateral Stick Inputs 


8 


Mach 0.75, Alt 20K - Right Stab Lock at 4 deg - Failure Inserted at t = 10 sec 





Figure 3: Aircraft pitch behavior along with the adaptive augmentation pitch acceleration signal for the right 

stab locked at 4 degrees at 10 seconds 
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Phi (deg) NN ad p (deg/s') 


Mach 0.75, Alt 20K - Right Stab Lock at 4 deg - Failure Inserted at t = 10 sec 





Figure 4: Aircraft roll behavior along with the adaptive augmentation roll acceleration signal for the right stab 

locked at 4 degrees at 10 seconds 
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NN ad r (deg/s ) 


Mach 0.75, Alt 20K - Right Stab Lock at 4 deg - Failure Inserted at t = 10 sec 





Figure 5: Aircraft yaw behavior along with the adaptive augmentation roll acceleration signal for the right stab 

locked at 4 degrees at 10 seconds 


11 


Mach 0.75, Alt 20K - Right Stab Lock at 4 deg - Failure Inserted at t = 10 sec 




Figure 6: Behavior of the performance errors (that drive the adaptation) along the three axes for the right stab 

locked at 4 degrees at 10 seconds 
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Can Pos (deg) Rud p os (deg) Ail Pos (deg) Stab Pos (deg) 


Mach 0.75, Alt 20K - Right Stab Lock at 4 deg - Failure Inserted at t = 10 sec 
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Figure 7: Aircraft surface commands axes for the right stab locked at 4 degrees at 10 seconds 
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Mach 0.75, Alt 20K - Canard Multiplier at -1 - Failure Inserted at t = 10 sec 



Figure 8: Aircraft pitch behavior along with the adaptive augmentation pitch acceleration signal for the canard 

multiplier set at -1 at 10 seconds 
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Mach 0.75, Alt 20K - Canard Multiplier at -1 - Failure Inserted at t = 10 sec 





Figure 9: Aircraft roll behavior along with the adaptive augmentation roll acceleration signal for the canard 

multiplier set at -1 at 10 seconds 
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Mach 0.75, Alt 20K - Canard Multiplier at -1 - Failure Inserted at t = 10 sec 



time (seconds) 
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Figure 10: Aircraft yaw behavior along with the adaptive augmentation roll acceleration signal for the canard 

multiplier set at -1 at 10 seconds 
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Mach 0.75, Alt 20K - Canard Multiplier at -1 - Failure Inserted at t = 10 sec 



-4 



Figure 1 1 : Behavior of the performance errors (that drive the adaptation) along the three axes for the canard 

multiplier set at -1 at 10 seconds 
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Mach 0.75, Alt20K- Canard Multiplier at -1 - Failure Inserted att= 10 sec 
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Figure 12: Aircraft surface commands axes for the canard multiplier set at -1 at 10 seconds 


V.Conclusions 

A stable discrete-time update law is investigated that provides augmentation to the existing controller only if the 
tracking error does not conform to its desired behavior. It is shown that not conforming to the desired first or second 
order behavior corresponds to the modeling error between the assumed and the actual plant. Thus the update law is 
modeling error driven, but in a direct adaptive framework. This update philosophy tries to regain the tuned 
performance that the control loop would have provided had there been no modeling error. This approach is being 
investigated for implementation as an alternate adaptive controller in the Intelligent Flight Control (IFC) architecture 
using the NASA Dryden high fidelity F-15 model. The results illustrate that we can achieve significant 
improvements in the aircraft behavior with the proposed approach. Future work will investigate a more 
comprehensive study of the approach including comparison of this approach for pilot handling qualities and frozen 
weight stability margin studies. We will also look at second order update methods for investigating the speed of 
convergence. 
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VI.Appendix 

Proof of stability for the update law given in Eg, (22) 

The weight error is defined as: 


W ad =w ad -w* 


Substituting in Eq. (25) gives: 


W ad .(£ + l) = W ad .(£) + 


a{k)*E j (k + 1)*P(£) 

'l + p r (W)] 


The performance error is given by Eq. (24) as: 

Efk + 1) = W^/pa) - W^/WP)*) 


The weight error update can thus be given as: 


w ad _a+i)=w ad ^)+ 


= ^a d; W - 
= ^a d; W - 


* [' w; di r p w - w ad . r wpw] * m 

) + p r (A:)P(A:)] 


a(k)* P(/t) 7 

' W ad ,(A:)-W; di ' 

*m 


i+p T (mk)_ 



a{k)*m T w aA jk)*m 

)+p r (£)pw] 


(Al) 


(A2) 


(A3) 


(A4) 
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The norm of the weight error is given as: 


W ad; (k + 1)|| 2 = w ad; (k + if W ad; (k + 1) 


(A5) 


Substituting Eq. (35) gives: 


,,2 I . a(k)P(kY W ad (Jt)p(Jfc) 

|w ad (^ + i) = w adi W- 


~l + p r (W)] 


W adi w- 


a(k)m T w ad imk)\ 

)+p r (£)p(£)] 


= W adj (£) y W adi (£) 


2a(k)w ddi (kf mm : 

r w dA ik) , 

* i 

m 2 m T m 

W adj (k) T $(k)$(k) T W ad . {k) 


\\+f{k)m_ 

T j 

)+p r (W)] 


\+tfm(k)_ 



= Iiw ad (rf +flW i w a ,(^p W p W -w a , W | 
-■ 11 i [i+p r (W)]J | [i+p r (W)_ 


(A6) 


W, 


ad, 


«(*)P(*) r P(*)l 

W adi (/c) r p(/c) 2 

“l + p r (^)P(^)]j 

\\+f{k)m_ 


(A7) 


For the condition given by Eq. (26), 


W ad .(A: + l) 2 - W ad (*)f < 0 


(AS) 


The weight vector error norm, being bounded and decreasing, converges to zero as k — > ®, This implies that 
W ad . ( k ) — > W adi _ Equation (38) also gives: 


lim 

k — >oo 


W ad . (k + 1) ” - W ad . (k) 2 ] = lim a(k) 

1 ' / &-» 00 


-2 + i 


a{k)${k) T ${k) 

1|| 

w ad x*fPW 

2 

- - 0 

~\+f{k)m] 

( 

\+f{k)m 

] 


w ad .W r P(*)|| 

' lim T V = 0 

*-*« l + p r (/t)P(/t)J 


(A9) 


lim E t (h) 1 = 0 

» oo 


Thus Ej (k) — > 0 as ^ oo 
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